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Abstract 

Colloidal particles are not simple rigid particles, in general an isolated particle is a system with many degrees of 
freedom in its own right, e.g., the counterions around a charged colloidal particle. The behaviour of model colloidal 
particles, with a simple phenomenological model to account for these degrees of freedom, is studied. It is found that 
the interaction between the particles is not pairwise additive. It is even possible that the interaction between a triplet 
of particles is attractive while the pair interaction is repulsive. When this is so the liquid phase is either stable only 
in a small region of the phase diagram or absent altogether. 

1 Introduction 

Theories for the behaviour of systems of particles usually apply to a model for the interaction between particles in which 
the energy of interaction is pairwise additive. The Hamiltonian is of the form 

N ' 

H(r N ) = ]T u(r«), (1) 

i,3'=l 

a sum of spherically-symmetric pairwise-additive potentials u{r). r N denotes the coordinates of the N particles and ry 
is the scalar separation of particles i and j. But experiments arc almost never on simple rigid particles. For example, the 
particles in a colloidal suspension are not simple rigid particles. Generally they are stabilised in the suspension either 
by being charged, in which case they are surrounded by a cloud of counterions, or by having short polymers grafted 
to their surface 0]. In either case we have not a rigid particle but a system with many degrees of freedom, e.g., the 
positions of the counterions, which fluctuate. Each particle is a system, with a free energy, susceptibilities etc., in its 
own right. Here we will treat systems of particles, each of which is a weakly fluctuating statistical mechanical system 
in its own right. Each particle will have a Landau-like free energy which couples to that of neighbouring particles. We 
will use the simplest possible form of this free energy, and show exactly and analytically that it leads to many-body 
attractions. These cannot be expressed as a Hamiltonian with the form of Eq. (0). 

This work is partly inspired by recent experiments on highly charged colloidal particles under conditions of minimal 
amounts of added salt. The potential of mean force between an isolated pair of colloidal particles with minimal salt 
concentrations has been measured and is purely repulsive [^). Yet the particles form crystallites which appear to be 
metastable at close to zero osmotic pressure || ; this is very hard to explain unless there is some sort of cohesive attraction 
between the particles in the crystallites. It does not seem possible to explain these findings with a Hamiltonian which 
depends only on the coordinates of the colloidal particles and has the conventional form of Eq. (|j) . Here we develop 
a simple phenomenological theory for colloidal particles which treats each particle as fluctuating statistical mechanical 
system. See Refs. § for recent theoretical work on understanding this behaviour starting from a Hamiltonian which 
explicitly includes the counterions and the electrostatic interactions. The theory is phenomenological as we simply 
assume that the state of a particle can be described by a single coarse-grained scalar variable, we do not derive this. 
This scalar variable fluctuates and these fluctuations are perturbed by the presence of another particle nearby. The 
perturbations due to the particles which surround any given particle add up, meaning that the state of a particle 
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changes with the number of its neighbours. Essentially, the more particles that surround any given particle the more the 
particle's fluctuations are biased towards values that minimise the interaction free energy and so the more attractive is 
the interaction between the given particle and all the surrounding particles. It is possible for two of our model particles 
to repel each other but for particles in clusters of more than two particles to attract each other. This can lead to 
coexistence between a dilute fluid and a dense crystal phase even when a pair of particles repel each other. We will show 
that generally if the fluctuations can be described by a scalar variable then the interaction between a pair of particles 
within a cluster of several particles is more attractive than between an isolated pair. 

In the next section we define our model. In section 3 we show that it exactly corresponds to a potential between 
structureless particle which contains both pair and triplet terms. Then in section 4 we apply a perturbation theory to 
obtain an approximation to the free energy of both the fluid and crystalline phases, which is accurate for long-ranged 
interactions between the fluctuating internal states. We show results for the phase behaviour and for the zero wavevector 
structure factor in section 5. Section 6 is a conclusion. 



2 Model 

Our model particles have an internal state specified by a single scalar variable. The value of the variable for a single 
isolated particle fluctuates weakly. Essentially, we view it as a coarse-grained variable |H obtained by averaging over 
some large number of degrees of freedom associated with each particle. An outline of how our mesoscopic Hamiltonian 
may be derived from a microscopic one is given in the Appendix. The interaction between a pair of the particles depends 
on the value of this variable. Thus, when two particles interact the mean values of the variables on the two particles will 
change; the particles 'polarise' each other, the interaction between them biases both internal variables towards values 
for which their interaction free energy is low. We included polarise in quotes because our model is phenomenological 
not electrostatic; cf Ref. || which considers the polarisation of one a charged colloidal particle by another. 

The effective Hamiltonian H of the system of particles has two parts: one part independent of the s variables, U, 
and the other is the free energy of the particles as a function of the TV internal variables, Fjy, 

H(r N ,s N ) = U(r N ) + F N (r N ,s N ), (2) 

where and s N symbolise the centre-of-mass and internal-variable coordinates, respectively, of all N particles. 

We need not specify U beyond saying that it should be such that the system has a well-defined thermodynamic limit 
|p^|] . Later on we will set U to be the sum of hard-sphere repulsions between pairs of the particles, 



N ' 

U(r N ) = -J2 (3) 

i,.7=l 



where U} ls is the hard sphere potential 



Uhs(r) = < (4) 

r is the separation of the centres of the interacting particles. The dash over the sum in Eq. (||) means that the sum 
is only over those terms for which j ^ i. Note that not all the double and triple sums below will have a dash, in the 
undashed ones terms with equal subscripts are summed over. 

Each particle has a dimensionless internal variable s. Now, we assume that the coupling between these variables on 
different particles is pairwise additive and spherically symmetric. Then the free energy Fn is, 



N N ' 

F N (r N ,s N ) = ^/W( Si ) + ^ J] f^i,*,^). (5) 

i—l 

is the free energy of a single isolated particle and is a function only of its s variable, is the difference in free 
energy between an isolated pair of particles and two isolated particles, and is a function only of the two s variables and 
the magnitude of the separation of the two particles. An isolated particle is a particle far from any other particle, an 
isolated pair of particles is a pair of particles far from any other particle. For weak fluctuations s is so we Taylor expand 
/(i) an d /(2) 

(s) = as 2 + a 3 s 3 + a 4 s 4 + • ■ • , (6) 
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where the linear term is missing as the variable s is defined so that if a particle is isolated its mean value is zero, and 

f (2 \s, s',r) = Mr) + Mr)(si + s') + ^ 2 (r)( S 2 + s' 2 ) + ^(r)ss' + • • • . (7) 
The coefficients a, etc. are derivatives of f^ 1 ' 

a= 2^^> (8) 

and «3 is 1/6 times the third derivative etc.. Similarly, the coefficients 4>o(r), $>i(?") etc. are derivatives of at a 
fixed separation of the particles. (f>o is the zeroth derivative, <pi is the first derivative with respect to either of the two s 
variables, 
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and similarly for the higher derivatives 4>2 etc.. 

Our Taylor expansions are quite general for particles which have a state which can be described by a single weakly 
fluctuating scalar variable and in which the interaction between two particles can be expressed in terms of this variable. 
If the variable is a vector or a tensor not a scalar then clearly there will be expressions analogous to that of Eq. (||) but 
vectorial or tensorial variables will in general lead to results very different from those found here. 

We truncate the Taylor expansions, Eqs. (||) and (Q), after their lowest nontrivial terms and substitute the resulting 
expansions into Eq. (^|) to obtain 



F N (r N ,s N ) = J2as 2 + \ £ [Mm) + Mm) (Si + Sj )\ . (10) 
1=1 i,j=l 

The inverse susceptibility, a, and the zeroth, (j>o, and first, <p%, order coefficients all have the dimensions of energy. 



So, for our Hamiltonian, Eqs. (|2j) and (10), to describe a system of particles accurately we require: (1) that the 
fluctuations in s of an isolated particle be sufficiently small that the higher order terms in may be neglected, (2) 
that the interactions between the particles can be described as a sum over interactions of pairs, (3) that this interaction is 
spherically symmetric, (4) that this interaction can be accurately described by function of a single scalar, coarse-grained, 
variable, and (5) that the interactions between fluctuations are sufficiently weak that the higher terms in /W and 
are small. 

For noninteracting particles the second sum of Eq. (|To|) may be neglected and Fn is just a sum of independent 
quadratic terms. The distribution of the s's is then Gaussian, which is correct for small fluctuations of a coarse-grained 
variable jl2|,|lj]. a is an inverse susceptibility of an isolated particle, the smaller it is the larger are the fluctuations. 
The interaction between a pair of particles is expressed as a Taylor expansion in the two s's truncated after the linear 
term. One particle feels an interaction due to another nearby particle, which couples to its order parameter s with a 



strength 4>i(r). The s-dependence of the free energy of Eq. (10) is what we would expect if the particles were weakly 
interacting, weakly fluctuating thermodynamic systems. 

3 Exact theory 

We start from the configurational integral Zn for N particles in a volume V and at a temperature T 

Z N = J dr N ds N exp{-H(r N ,s N )/kT), (11) 

with the Hamiltonian given by Eqs. (§) and (^p|). The Helmholtz free energy A is then 

A = -kT\n(Z N A N /M) , (12) 
where A derives from the integration over the momenta. Using Eqs. (||) and ( |To| ) in Eq. (|Tl|), 



J dr N exp(-U(r N )/kT) J ds N exp ^- J^(a/kT)s 2 - ± [Mm)/ kT + <f>i(m) ( s i + s j) / kT ]j 



^ i-^-i £ ^) n {/>- -t^-\>. <«) 

i,j=i j i=i ^ y j=i jjti 
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where to obtain the second line we expressed the integral over the s variables as a product of integrations over each one 
and then grouped all the terms which depend on each Si together. Each integration over an s variable is independent 
and can be done easily 
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(14) 



where to obtain the second from the first line we expressed the square of the sum as a double sum. To obtain the third 
line we converted the product of exponentials to the exponential of a sum. The factor in front of the integration is of 
course very familiar: it is just the partition function of N independent simple harmonic oscillators. It does not depend 
on density and so has no effect on the phase behaviour. Note that in the triple sum although neither j nor k can be 
equal to z, j can be equal to k. The restrictions on j and k derive from the fact that a particle cannot interact with 
itself, which would correspond to j,k = i but as both j and k in the triple sum come from the square of a single sum 
they are in effect from the same interaction and therefore are allowed to be equal. We can rewrite Eq. ( |l4| ) by extracting 
the j = k terms from the triple sum, 
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JV 



AakT 
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where the triple sum has a dash to indicate that terms in which any of the three subscripts are the same are excluded. 
The configurational integral, Eq. ([l5|), is an integral only over the positions of the N particles. Neglecting the irrelevant 
(for the phase behaviour) prefactor it is nothing but the configurational integral of the Hamiltonian, H e ff, 



N 



Heff(r N ) 



U(r N ) + 5 E 



Mr*)* 



2a 



JV 



i,j,fe=l 



[<fti(T-jj)(fti(rifc) + <h^Tij)^£jk) + <f>i(nk)<fii(rjk)} 
2a 



(16) 



where we have rewritten the summand of the triple sum to make it symmetric with respect to the three indices. The 
phase behaviour of our model particles will be identical to that of structureless, spherically symmetric particles with 
interactions described by the Hamiltonian H e f j of Eq. (|l^) . The Hamiltonian is that of a triplet or 3-body potential; 
the first two sums are conventional sums over pair potentials but the last sum is over a 3-body potential. We started 
with a pair potential which depended on the s variables as well as the positions, integration over the s variables resulted 
in a potential which depends only on the positions but is no longer a simple pair potential. 

The effective Hamiltonian, H e ff, can lead to behaviour qualitatively different from that of a Hamiltonian which is a 
simple sum over a pair potential. To see this we will compare the interaction between a pair of particles to that between 
a triplet of particles. For a pair of particles, 1 and 2, a distance r apart, the interaction is from Eq. (|l6|), 



ff o// (H) = tf(H) + to(r)-- 



iMr) 2 



(17) 



Note that the contribution due to the fluctuations, the part inversely proportional to a is always negative regardless 
of the sign of tf>i, fluctuations always produce an effective attraction. Now, consider 3 particles, at the corners of an 
equilateral triangle of side r for simplicity. We denote the set of 3 coordinates marking the corners of an equilateral 
triangle by r^. The interaction is, from Eq. (|l6|), 



H eff (T 3 e ) = U(4)+3Mr)-3 



Mr) 2 



a 



(18) 
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If we compare this with the interaction of a pair, Eq. (fi~7|), we see that the ratio of the 4>\/a to (f>o term has doubled. 
For a pairwise additive potential, the interaction free energy would be simply three times Eq. ([l7|): 3</>o — {Z/2)<j>\/a. 
The relative contribution from the fluctuations in the s variables has doubled. This contribution is always attractive 
so the free energy of attraction is always more negative than for a pairwise additive potential. This is a general result 
for a weakly fluctuating scalar variable. Our result for the interaction between three particles may be compared with 
the more microscopic work on the interaction between three particles of Lowen and coworkers. They considered triplet 
interactions between charged colloidal particles jL5) and between star polymers [jl6| . 

When (f>o{r) > and </>o(r) < <j>i{r) 2 /a < 2(j) (r), the interaction minus the part from U is repulsive, i.e., greater 
then zero, for a pair r apart but attractive for a triplet at the corners of an equilateral triangle of side r. If U is 
some short-range repulsion, e.g., a hard sphere repulsion, which is zero at some sufficiently large value of r, then it is 
possible for a pair of our model colloidal particles to repel each other, but for a triplet of them to attract. For other 
arrangements of 3 particles the sign of the interaction free energy will depend on the values of the 3 separations of 
the centres of the particles and on the distance dependence of 4>o and <j)\. However, this should not obscure the basic 
fact that the interaction free energy of 3 particles can be negative when the interaction of 2 particles is positive. For 4 
particles at the corners of a tetrahedron the interaction minus the part from U is 6(/>o — 901/ a Again, it is possible to 
obtain negative interaction free energies even when the interaction free energy of a pair is always positive. Indeed it is 
possible to obtain negative free energies even when they are positive for three particles at the corners of an equilateral 
triangle. 



4 Mean-field free energy 

In order to able to keep the theory simple we now specialise to an interaction between the internal variables which is both 
weak and long-ranged. We also set U to be the hard-sphere interaction. Weak in the sense that |0 o ( r )l: \4>i{ r )4>i( r ') / a \ ^ 
kT for all values of r, r' , and long-ranged as both 4>o(r) and <j)\{r) decay to zero over a characteristic length scale which 
is much longer than the hard-sphere diameter a. So the interaction Hamiltonian, Eq. (fl6|), which is a function of the 
positions only, consists of a hard-sphere interaction, U, and a weak long-ranged interaction which has two parts: a pair 
potential and a triplet potential. Weak but long-range many-body interactions have been considered by the author in 
Ref. @. 

We view the long-range part of the interactions as a perturbation 111] . The free energy is expressed as that of hard 
spheres, A} lSl plus the difference in free energy between a system of our particles and that of hard spheres, AA, 

A = A hs + AA. (19) 

For the hard sphere free energy we use the expression of Carnahan and Starling for the fluid phase and the fit 
to simulation data of Hall [^(J for the solid phase. We approximate AA by an average of the perturbing part of the 
Hamiltonian, H e jf — U, 

AA = I dr N (H eff (r N ) ~ U(r N )) exp (-U(r N )/kT) 

Zhs 

where Z^s is the configurational integral for hard spheres. 
Substituting Eq. dl6| ) into Eq. (po|), 

AA = Z~ h l J dr- exp (-U(r»)/kT) [l N (N - 1) ^ (r 12 ) - M^) - - 1)(N - 2) gM^iM j . (21) 
The n-particle density of hard spheres, p^J , and their n-particle distribution function, gj^J , are defined as Jl8| 

N\ Jexp[-pU(r N )]dr N - n 



(22) 



(N-n)\ Z hs 

The 1-particle density, (r) , is not assumed to be uniform so that the theory applies to crystalline as well as fluid 
phases. Using, Eq. (p2[) in Eq. (|2l|), we obtain 



AA = I J &v*p%{v u r 2 ) Uo(r 12 ) " ^ ^(^2^3)^2)^13). (23) 
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As <f>o,4>i decay to zero only when the separations of the particles are much larger than cr, the integrals in Eq. (|23|) 
are dominated by configurations when the spheres are far apart, i.e., where the pair separations are much larger than 
a. In a fluid the one particle density is a constant, p^ = p, and at separations large with respect to a the distribution 
function is close to one, g^J ~ 1. Thus in the fluid phase the integrands of Eq. ( |23| ) can be simply approximated by 
P 2 {4>o( r i2) — 4>i( r i2) 2 /2a;) and p 3 4>\{r 12) 4>i(tiz)- In a crystalline phase, although there are long-range correlations in 
Phs t ne one particle density, pj, s , averaged over a unit cell is just p. The attractive interaction between particles has a 
range much larger than the lattice constant of the lattice (a little larger than a) and so as 4>q and <j}\ vary little across a 
unit cell we can regard (r) as approximately constant at its average value, p. Similarly, as we change any one of the 

(n) 

n position vectors upon which the n-body distribution function, p hs , depends the density oscillates rapidly over each 

unit cell but averages to p. So, in the crystalline phase as well as in the fluid phase we approximate p^ by p n . Then 
the integrands are the same as in a fluid phase. We have 

2 \ i / poo \ 2 



AA = -Vp* J dr ^ (r) - I - — V p A (J dr^(r) j . (24) 

Now the term coming from the square of <f>\ is negligible if 4>i is long-ranged. To see this consider an explicit functional 
form for <j>\. We choose a Kac potential, 

(j>i(r) = e7 3 exp(— jr/a), (25) 

The range of this function is cr/j, i.e., j^ 1 is a range in units of a. Now the integral of Eq. (|25| ) over 3-dimensional 
space is 0(ea 3 ); essentially the integrand is of order ey 3 within a volume of order (cr/7) 3 , and negligible outside of this 
volume. However, the integral of the square of Eq. (|25| ) over 3-dimensional space is 0(ecr 3 7 3 ). For a long-range (j>i, the 
inverse range 7^1 and so the integral of the square is negligible. We chose a specific functional form for 4>i simply for 
clarity, our conclusion that the integral of the square is negligible holds for any long-range slowly decaying function. 



So, in Eq. (24), we neglect the <j)\ term, integrate over the remaining terms in Eq. (E4I) and obtain 



AA 1 \{p Vl f 

~W = 2 ~ 4 a ' (26) 

where and z/i are the integrals 

roc 

dr^(r) =v h t = 0,l. (27) 



The pressure p can be easily derived from the free energy, Eq. (^6|) , 

1 , 1 p 3 v} . . 

P = Phs + t;P "o~t; — , 28 

2 la. 

where phs is the pressure of hard spheres. The contribution from the fluctuations in the s variables to the pressure is 
cubic in the density as it must, it is equivalent to a long-range 3-body attraction ]T^ ]. The chemical potential p is then 
simply given by p = A/N + p/ p. The pressure and chemical potential as functions of the temperature and density allow 
us to calculate phase diagrams. 



5 Results 



The free energy, Eq. (26), depends on a single parameter, the dimensionless ratio R = v^aa jv\. It ranges from — 00 
to +00; the sign of R is determined by the sign of Pq as a must be positive. The larger the magnitude of R the more 
dominant is the part of the interaction potential which does not depend on s and so the closer the interaction is to 
a pairwise additive potential. In the R — > — 00 limit the interaction is a simple pairwise additive attractive potential. 
Then the free energy is just a modification of the free energy of van der Waals fif , ^l|-p3|. The modification being the 
replacement of his approximate free energy of hard spheres by that of Carnahan and Starling fl9[ | in the fluid, and that 
of Hall [E0| in the crystal. In the R — » +00 limit the potential is again pairwise additive but it is repulsive. A long-range 



repulsion tends to cause not bulk vapour- liquid separation but microphase separation, see Refs. |22, 24 1. The R — > ±00 
limits result from either the interaction between particles not depending on the value of the s variables, vi — > 0, or the 
fluctuations of the s variables tending to zero, a — > 00. 
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When R is finite then the interactions depend on the fluctuations in s and are no longer pairwise additive. For R < 
the interactions are attractive even in the absence of fluctuations in s; the fluctuations in s merely make the attractions 
stronger and not pairwise additive. For R > the pairwise interactions are repulsive while the 3-body interactions 
from the fluctuations are, as always, attractive. The part of the free energy due to fluctuations, the last term in Eq. 
(pq), is one power higher in the density than the part independent of the fluctuations, the second in Eq. (p6|). Thus, if 
vq > the interactions are always repulsive at sufficiently low densities but become attractive, i.e., contribute a negative 
amount to the free energy of the system, at a density 2R/<r 3 . 

We now discuss some example phase diagrams. The diagrams are in the density-temperature plane. We use the 
reduced density r\ — (7r/6)per 3 , which is equal to the fraction of the volume occupied by the hard cores, and a reduced 
temperature T* which is either kTa- 3 /\vo\ or kTaa 6 jv\. For reference we plot the familiar phase diagram of hard 
spheres with a long-range attraction, often called the van-der-Waals fluid |^l|-|23|]; the R — > — oo limit. This has a free 
energy given by Eqs. (|l^) and ((2^) with v\ = and vq < 0. The first condition means that the internal variables on 
different particles do not couple and so the attraction is a simple pair potential and the second condition makes the 
long-range interaction attractive so there is a vapour- liquid transition, i.e., phase separation into two fluid phases of 
different densities. The phase diagram is shown in Fig. [l; there is a large temperature range over which there is stable 
vapour-liquid coexistence. This is of course not new, in the v\ = limit our model reduces to the most basic model of 
particles which form a liquid. 

In Fig. H we show the phase diagram of particles with R = 0. The long-range interaction is purely proportional to 
the internal variable s, i.e., z/ = 0. There is vapour-liquid coexistence over a range of temperatures but this range is 
much smaller than for the simple pair potential of Fig. |l|. Also, the density at the critical point is higher; it is almost 
double that in Fig. [j]. The interaction between s variables is effectively a long-range 3-body attraction, i.e., an attraction 
which is van-der-Waals like except for the fact that it is between triplets not pair of particles . The phase diagram, 
Fig. |[ differs from Fig. 2 of Ref. jl7j only in that the temperature scale is different. 

In Fig. H we show the phase diagram of particles with R — 2. There is no vapour-liquid coexistence at equilibrium, 
only a fluid-crystal coexistence region which broadens dramatically at low temperatures. The liquid phase has disap- 
peared. For other examples of liquid phases disappearing see Ref. [p5| ; the most studied system in which the liquid 
disappears is that of particles with a short-range, pairwise additive attraction |^6|-^]. Within the fluid-crystal coexis- 
tence region the pressure and chemical potential have van-der-Waals loops so our bulk free energy predicts vapour-liquid 
coexistence within the fluid-crystal coexistence region. Note that the density at the critical point is even higher than in 
Fig. | 

We have only considered fluid phases and crystalline phases. However, sufficiently strong long-range repulsions 
can transform a bulk phase separation into microphase separation j2^, ^4j. Essentially, microphase separation occurs 
when the interaction between particles is repulsive at the largest separations at which they interact separately and the 
repulsion is sufficiently strong and long ranged. Thus our present theory, which neglects the possibility of microphase 
separation, should not be used if the long-range interaction is predominantly repulsive, i.e., if vo > 0, and either </>o is 
longer ranged than (f>i, or R S> 1. 

5.1 The structure of the fluid phase 

We have shown that the phase behaviour of our particles can be very different from that of a simple pair potential but 
what about the structure? Of course the pair distribution function or equivalently the structure factor | fl8[ depends on 
the functional form of the interactions, i.e., on the precise forms of the functions </>o( r ) and 4>\{r). However, we have 
only considered long-range interactions and there the bulk phase behaviour is insensitive to the precise details of <po ( r ) 
and 4>i(r), it only depends on their integrals vq and i>i. Also, long-range interactions only affect the structure factor 
S(q) at small wavevectors q; small meaning of order the reciprocal of their range or less. Thus, we will consider only the 
zero wavevector limit, 5(0), of the structure factor. This is simply related to a thermodynamic quantity, the isothermal 
compressibility xt, by [[L8| 

S(0)=pkT XT , (29) 



where 



Using Eq. (08) for the pressure 



tf'-'d),- (30) 

(31) 
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Eqs. (29) and (|3lj) give us the zero wavevector structure factor of our particles. Below we will consider the density 
dependence of 5(0). It will therefore be useful to obtain 5(0) as a density expansion. To do this we start by inserting 
the virial expansion for the pressure into the definition of the isothermal compressibility xt, Eq. (|30"|). Inserting this 
expansion into Eq. ( p9| ) for 5(0) yields 

5(0) = 1-2B 2P + 0( P 2 ), (32) 
the initial slope of 5(0) is equal to minus twice the second virial coefficient B^. For our particles 

B 2 = ^a 3 - ^o, (33) 

where the first term on the right hand side is the second virial coefficient of hard spheres. Note that the fluctuations in 
the s variables do not contribute to the second virial coefficient. This is only the case in the long-range interaction limit. 
In section 3 and in Eq. ( pi| ) we found that in general the fluctuations did contribute to the interaction between pairs 
and hence to the second virial coefficient but that this contribution was weaker than its contribution to the interaction 
between larger numbers of particles and to the third virial coefficient. The term in the free energy per unit volume 
proportional to v\ /a which comes from the fluctuations varies with density as p 3 and so it contributes (only) to the 
third virial coefficient. 

In Fig. ^ we have plotted the zero wavevector structure factor for a van-der-Waals fluid, hard spheres plus a 
long-range pairwise additive attraction (the solid curve), and for an attraction which is proportional to s (the dashed 
curve). They are both at temperatures just above the critical temperature. The peaks in 5(0) are due to the nearby 
critical points, where 5(0) diverges of course. The striking difference between the curves is their opposite slopes at low 
densities. The limiting slope at vanishing densities is equal to — 2i?2, Eq. (|32|). For a van-der-Waals fluid the second 
virial coefficient is negative at the critical temperature so the slope of 5(0) is positive until above the critical density. 
However, interactions between fluctuations do not contribute to the second virial coefficient and so it is equal to that 
of hard spheres which is positive. Thus the slope of 5(0) is negative at low densities and 5(0) goes through a minimum 
below the critical density. 



6 Conclusion 

Particles which have a fluctuating internal state interact in a way which is different from structureless particles. The 
interaction is intrinsically many body and cannot be described using a pair potential. It is even possible for pairs of 
particles with a fluctuating state to repel but larger clusters of particles to attract each other. We have studied a rather 
generic model of a particle with a fluctuating state. Remarkably, because of the very simple form of our Hamiltonian's 
dependence on the variables which describe the states of the particles, see Eq. (|Io|), we have been able to integrate 
over these variables exactly and analytically. The result is an interaction between the particles which has both pair 



and triplet interactions, Eq. (16). Our model only included the leading order terms in Taylor expansions for the free 
energy associated with the internal states of the particles. It is possible to include higher order terms although then 
the s variables will have to integrated over approximately or numerically, alternatively computer simulation could be 
used. In any case the effective potential will then contain not only pair and triplet terms but terms of all orders, 4-body, 
5-body, etc.. 

It is perhaps of interest to contrast the particles with mesoscopic fluctuations studied here to the rods with microscopic 
fluctuations studied by Ha and Liu see also Refs. ]3l|-|3l| . They studied the interactions between parallel rods and 
allowed the charge density along each rod to fluctuate; if a rod is along the z-axis then we can define a 1-dimensional 
charge density p(z), which will fluctuate. The Coulomb interaction between charges will then couple the fluctuations 
in p(z) in nearby rods. Ha and Liu found that not only can this coupling generate an attraction but that it is strongly 
nonpairwise additive. In contrast with their model our is rather general and simpler. The results for our rather non- 
specific model show clearly that these findings of Ha and Liu are rather generic. The disadvantage of our model with 
respect to that of Ha and Liu is that to relate it to experiment requires determining or at least guessing the appropriate 
values of the parameters of our phcnomcnological interaction. 

We hope that our model will be useful for determining generic differences in behaviour between rigid particles 
interacting via a pair potential and particles with a fluctuating state. For example, although all our calculations have 
been for a single component they are easily generalised to mixtures. When this is done, it is straightforward to see that, 
if the fluctuations between the internal states of particles of different components are only weakly coupled then this will 
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tend to drive phase separation of these components. However, differences between the inverse susceptibility, a, cannot 
drive phase separation. At least within our theory which includes only leading order terms. 

It is a pleasure to acknowledge discussions with B.-Y. Ha, Y. Levin, R. van Roij and P. B. Warren. 



Appendix 

Rather generally, when we are evaluating a configurational integral, with the aim of calculating the free energy, we 
can split the variables over which we are integrating into two sets. Then integrating over one set is always possible in 
principle, and it results in an effective Hamiltonian which is a function only of the remaining set of variables. However, 
this effective Hamiltonian is not in general pairwise additive, simply because there is no reason for it to be. If we denote 
the set of all variables by v N+M , and the two subsets by v N and v M then the configurational integral is 



Z = J dv N+M cxp (-H (v N+M )/kT) = J dv N exp (~H'(v N )/kT) , (34) 

where 

exp (-H'(v N )/kT) = J dv M cxp (-H(v N+M )/kT) . (35) 

H' will not in general be expressible as a sum over a pair potential even if H can be. Physically, this is most often useful 
when the two sets of coordinates, v N and v M , are very different. We consider a model in which the set of coordinates 
v M is coarse grained |12]] , i.e., the integration over v M is left as a function of a set of coarse-grained variables, so instead 
of Eq. ([35]) we have, 

cxp (-H'(v N , c L ))/kT) = [ dv M exp (-H(v N+M )/kT) JJ S (c a (v M ) - c a ) , (36) 

•* a=l 

where c L is a set of L <C M coarse-grained variables. c a is the ath coarse-grained variable and c a (v M ) is the definition 
of the ath coarse-grained variable in terms of the microscopic coordinates, v . The configurational integral is then 

Z = J dv N dc L exp (-H'{v N ,c L )/kT) . (37) 



Our model has a configurational integral, Eq. ([II]) , of this form, with L = N as there is one coarse-grained variable per 
particle. The v N and c N coordinates are the positions and s variables of the particles. 
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Figure 1: The phase diagram of hard spheres plus a long-range pairwise additive attraction, a van-der-Waals fluid. 
v\ =0, vq < and the reduced temperature T* — kTcr 3 / \uq\- The thick solid curves separate the one and two-phase 
regions. The letters V, L and C denote the regions of the phase space occupied by the vapour, liquid and crystalline 
phases, respectively. The horizontal thin lines are tie lines connecting coexisting densities, and the dotted line connects 
the three coexisting densities at the triple point. 
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Figure 2: The phase diagram of hard spheres plus a long-range interaction proportional to the fluctuating internal 
variable s. i/o = 0, and the reduced temperature T* = kTaa 6 jv\. See caption to Fig. |l| for the meaning of the letters, 
curves and horizontal lines. 
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Figure 3: The phase diagram of hard spheres plus a long-range repulsion and a long-range attraction proportional 
to the fluctuating internal variable s. R = —0.5, and the reduced temperature T* = kTaa e jv\. The letters F and C 
denote the regions of the phase space occupied by the fluid and crystalline phases, respectively. The dotted curve is the 
bulk free energies prediction of a vapour-liquid-like transition within the fluid-crystal coexistence region. See caption 
to Fig. |j] for the meaning of the other curves and the horizontal lines. 




Figure 4: The zero wavevector structure factor of our particles in the fluid phase, >S'(0), as a function of the volume 
fraction r\. The solid curve is for an attraction which does not depend on s; the potential with the phase diagram 
Fig. [lj Po/(a 3 kT) = —10, which is just above the critical temperature The dashed curve is for an attraction which is 
proportional to the fluctuating internal variable s; the potential with the phase diagram Fig. ^. v\j(a^kTa) = 19.4, 
which is again just above the critical temperature 
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